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We studied underdoped high Tc superconductors using a spinon-dopon approach (or doped-carrier 
approach) to t-t'-t"-J model, where spinon carries spin and dopon carries both spin and charge. 
In this approach, the mixing of spinon and dopon describes superconductivity. We found that 
a nonuniform mixing in fc-space is most effective in lowering the t' and t" hopping energy. We 
showed that at mean-field level, the mixing is proportional to quasiparticle spectral weight Z-. 
We also found a simple monte-carlo algorithm to calculate Z- from the projected spinon-dopon 
wavefunction, which confirms the mean-field result. Thus the non-uniform mixing caused by t' and 
t" explains the different electron spectral weights near the nodal and anti-nodal points {i.e. the 
dichotomy) observed in underdoped high Tc superconductors. For hole-doped sample, we found 
that Z is enhanced in the nodal region and suppressed in the anti-nodal region. For electron doped 
sample, the same approach leads to a suppressed Z in the nodal region and enhanced in the anti- 
nodal region, in agreement with experimental observations. 

PACS numbers: 71.10.-w, 74.72.-h, 74.25. Jb 



I. INTRODUCTION 

One powerful experimental technique to study high-Tc 
material is the Angular Resolved Photoemission Spec- 
troscopy (ARPES)[1]. ARPES study for the pseudogap 
region showed a strong anisotropy of the electron spec- 
tral function in momentum space [2, 3]. Basically it was 
found that in the nodal direction, excitations are more 
quasi-particle like; while in the anti-nodal direction, exci- 
tations have no quasi-particle peak. This is the so-called 
dichotomy. If one lowers the temperature to let the ma- 
terial to go into superconducting phase, it was found that 
anti-nodal direction also has a small quasi-particle peak. 
Tunneling experiments show that the underdoped sam- 
ples are inhomogeneous[4-6]. Due to this inhomogeneity, 
it is possible that the underdoped sample can be sep- 
arated into optimal doped regions and underdoped re- 
gions, and the quasi-particle peak only comes from the 
optimal doped region. With such a point of view, it is 
possible that even the superconducting phase can have a 
very anisotropic electron spectral function in momentum 
space. 

Exact diagonalization on t-t'-t"-J model {t' and t" 
stand for next nearest neighbor and next next nearest 
neighbor, respectively) with 32 sites has been done[7, 8] 
for hole-doped case (one hole doped). It was found that 
if t' = t" = 0, J = 0.3i, then the quasiparticle weight 
Z_ is almost a constant along the direction (tt, 0)-(0, tt): 
Z_ = 0.311 at (7r/2,7r/2), and Z_ = 0.342 at (7r,0). 
However if one put in t' = -0.3t,t" = 0.2i, J = 0.3t, 
which is an optimal parameter fitting for Sr2Cu02Cl2, 
then there is a strong dichotomy feature: Z_ — 0.353 at 
(7r/2,7r/2), and Z_ = 0.029 at (7r,0). This suggests that 
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the dichotomy can be a result of t' and t" hopping. 

Exact diagonalization was also done for the electron 
doped case (a few electrons doped on 32 sites) [8], where 
was measured. Due to the particle-hole symmetry at 
half- filling, we know that if t' = t" = 0, Zj^. and Z_ are 
equal up to a momentum shift of (tt, tt). Therefore Z^ 
are also flat along the direction (tt, 0)-(0, tt) in pure t-J 
model. But when we put in t' = —0.3t,t" = 0.2t,J = 
0.3i, the particle-hole symmetry was broken. was 
found to develop a strong anisotropy along (tt, 0)-(0, tt): 
Z+ = 0.005 at (7r/2,7r/2), and Z+ = 0.636 at (7r,0). 

What mechanism can destroy the quasi-particle co- 
herence in the anti-nodal region? The simplest thing 
comes into one's mind is that we need some other things 
to destroy it. For example, neutron scattering experi- 
ments indicate that there are some low energy magnetic 
fluctuations[9-12], and it was proposed[2, 3] that magnon 
scattering process can destroy the quasi-particle coher- 
ence in anti-nodal region. In this paper, however, we 
propose a physically different scenario: The dichotomy 
is due to the t' and t" hopping terms. The quasiparti- 
cle spectral weight Zk is naturally suppressed in some 
region in fc-space to lower the t' and t" hopping energy. 
This contradicts a naive thinking that hopping always 
enhance Zk. Using the t-t'-t"-J model, we will show that 
the new scenario can explain the distribution of Zk for 
both hole-doped and electron-doped samples in a unified 
way. 

If we believe that the dichotomy is driven by the t' 
and t" hopping terms, then there is an important issue: 
Is there a mean-field theory and the corresponding trial 
wavefunction that captures this mechanism? 

One way to understand high-Tc superconductors is to 
view them as doped Mott-insulators. Under Zhang- Rice 
singlet mapping[13], the minimal model which includes 
the essential Mott physics is t-J model on square lattice. 
On the analytical side, a powerful mean-field theory for 
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t-J model, the slave-boson approach, was developed [ 14- 
16]. This approach emphasizes the fractionalization pic- 
ture of the doped Mott insulator: electron is splited into 
a spinon (a fermion with spin and no charge) and a holon 
(a boson with charge and no spin), which characterize the 
low energy excitations of the doped Mott insulator. This 
mean-field approach also successfully predicted the pseu- 
dogap metal for underdoped samples. On the numerical 
side, the same physics picture gives rise to the projected 
BCS wavefunction[17](pBCSwf), which turns out to be 
a very good trial wavefunction for t-J model. However, 
more detailed studies of pBCSwf[18, 19] indicate that the 
slave-boson approach fail to explain the dichotomy. So 
a momentum dependent quasiparticle weight remains 
to be a big challenge for slave-boson theory. 

In this paper, we will use a new spinon-dopon 
approach [20] and the corresponding trial wavefunction 
to study the underdoped samples. Instead of using 
spinous and holons, in the new approach, we use the 
spinous and the bond states of spinous and holons to de- 
scribe the low energy excitations. The bond states of 
spinous and holons are called dopons which are charge-e 
spin-i fermions. The spinon-dopon approach leads to a 
new trial wavefunction, the projected spinon-dopon wave 
function (pSDwf). The new trial wave function turns out 
to be an improvement over the old projected BCS wave- 
function (pBCSwf). 

The holon condensation in slave-boson approach corre- 
spond to spinon-dopon mixing. However, in the spinon- 
dopon approach, the mixing can have a momentum de- 
pendence, which is beyond the mean-field slave-boson 
approach. If we set the mixing to have no momentum 
dependence, then the pSDwf turns out to be identical to 
the old pBCSwf. So the pSDwf is a generalization of the 
pBCSwf. 

Now the question is, why the mixing wants to have 
strong momcntiim dependc^nce? The answer is that 
the wavefunction with momentum dependent mixing can 
make the hopping more coherent, and therefore gain hop- 
ping energy. Roughly speaking, the pSDwf with momen- 
tum dependent mixing is the summation of the old pro- 
jected BCS wavefunction together with hopping terms 
cjcj acting on it. Here one should notice that the old 
pBCSwf, with uniform mixing, already has a pretty good 
t hopping energy. But to have a good and hop- 
ping energy, the mixing needs to have a momentTim de- 
pendence, along the direction from (tt, 0) to (0,7r). Our 
Monte Carlo calculation shows that pSDwf with momen- 
tum dependent mixing is indeed a better trial wavefunc- 
tion in energetic sense. To get a quantitative sense how 
big is the improvement, we find that the energy of a 
doped hole in pSDwf is about OAt lower than that of a 
doped hole in pBCSwf. This is a very big improvement, 
indicating that the spin-charge correlation (or more pre- 
cisely, the spin configuration near a doped hole) is much 
better described by pSDwf than pBCSwf. 

Can one measure the momentum dependence of mix- 
ing? The answer is yes. In this article we will show that 



the mixing is directly related to the quasi-particle 
weight, which is measurable in ARPES. Roughly speak- 
ing, mixing is proportional to We have also devel- 
oped a Monte Carlo technique to calculate Z- . The cal- 
culation shows that momentum dependent mixing pSDwf 
indeed has strong anisotropy in momentum space, and 
consistent with the observed dichotomy. 

Comparing pBCSwf and pSDwf, we like to point out 
two wave functions have similar background spin-spin 
correlation and similar spin energy. However, that pBC- 
Swf does not capture the detailed charge dynamics. 
The new trial wavefunction, pSDwf, contains more cor- 
rect spin-charge correlation. As a result, the energy 
of doped holes/electrons is much lower in the pSDwf. 
The holes/electrons in the pSDwf reproduce the correct 
momentum dependence of quasi-particle spectral weight. 
We also expect our pSDwf to have a strong momentum 
dependence in quasi-particle current, which may explain 
the temperature dependence of superfluid density of High 
Tc superconductors [21]. 

II. SPINON-DOPON APPROACH AND PSDWF 

A. Slave-boson Approach and Projected BCS 
Wavefunction — Why the approach fails to capture 
A;-dependent features? 

Why would we want to introduce the spinon-dopon 
approach to t-J model? Let us firstly look into the pre- 
vious mean field approach, more specifically, slave-boson 
approach. The general t-J model can be written in terms 
of electron operator: 

Htj =J ^ iSi- Sj - ^UiTij j 

-^UjP{4c,+c]ci)p. (1) 

ij 

Here the projection operator P is to ensure the Hamil- 
tonian is acting within the physical Hilbert space: one 
each site, the physical states are ] t), I i) or ]0), i.e., no 
double occupancy. 

Slave-boson approach[14, 15] emphasizes the spin- 
charge separation picture. In that approach, one splits 
electron operator into spinon and holon operators: 

Cia=fiabl (2) 

where / is spinon, carrying spin 1/2 and charge 0, i labels 
site and a labels spin; b is holon, carrying spin and 
charge 1. This splitting enlarges the Hilbert space. To 
go back to physical Hilbert space, a local constraint is 
needed: 

fyn+fyii+bibi=i. (3) 

Due to spin interaction, spinous form a d-wave paired 
state. The superconducting phase is realized through 
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an additional holon condensation at momentum fc ~ 0. 
Within such a construction, the quasi-particle weight Z is 
proportional to doping x everywhere in k-space, in both 
nodal and anti-nodal region. To see this, one can simply 
look at the mean-field Green function of electron: 



(cfe4)-(Co&fe=oA.At> = -^(.44^> 



(4) 



Therefore x is the residue of quasi-particle pole and Z = 
X is independent of k. 

Slave-boson approach is supposed to capture the 
physics of spin-charge separation. It has successfully 
generated the phase diagram of High-Tg superconductor. 
But this approach, at least at mean-field level, could not 
capture some more detailed features, such as momentum 
dependence of quasi-particle weight or the quasi-particle 
current. One can argue that including gauge fluctuation, 
those detailed features may be reproduced, but here we 
will try to develop another approach which can capture 
these features at mean-field level. 

Before we go into the new spinon-dopon approach, let 
us see how far one can go using slave-boson approach. 
One can actually try to build a trial wavefunction based 
on slave-boson mean-field approach. We know that the 
mean-field approach enlarged the Hilbert space, and the 
resulting wavefunction lies outside the physical Hilbert 
space. Only when one includes the full gauge fluctuations 
can one go back to the physical Hilbert space. 

So one way to include the full gauge fluctuations, is 
to build the mean-field ground state first, and then do a 
projection from the enlarged Hilbert space to the physical 
Hilbert space. The wavefunction after projection would 
serve as a trial wavefunction for the physical Hamilto- 
nian. This projected wavefunction is supposed to incor- 
porate the effect of gauge fluctuation of the slave-boson 
approach, and may answer the question that, after in- 
cluding gauge fluctuation, whether slave-boson approach 
can capture the detailed features like dichotomy. 

The mean-field ground state for underdoped case can 
be constructed as follows. Let be the number of holes, 
N f is the number of spinous, and N = Nh + Nj is the 
total number of sites. The slave-boson mean-field ground 
state is then given by 

|$SB,„> = {bU,f- Hiu, + Vkfl^fi^^m, (5) 
k 

where the spinon part of the wavefunction is a standard 
d-wave pairing state: 



Vk 



Aik) 



(6) 



where 



= -2x(cos/c,,: + con ky) - II 
A{k) = A{cos kx — COS ky) (d-wave). (7) 

Here /i is the chemical potential to give the correct av- 
erage number of spinon fi^fia) = ^f] X ^tnd A 



are mean-field parameters which have been found to be 
Y = 2 [22] at half-filling, and ^ decreases to zero at 
doping around J/t. 

Now one can do a projection to go back to the physical 
Hilbert space. The constraint for physical Hilbert space 
is Eq.(3). This constraint ensures that the total number 
of spinon must be Nf and there is no double occupancy 
of spinon: spinon number Uf^t at site i has to be cither 
or 1. One can easily see that the resulting wavefunction 
is the usual Projected d-wave BCS Wavefunction (pBC- 
Swf): 



Uk 



V]^Ci,t^ c 



k^'--ki)\^) 

Nf/2 



|0), 



(8) 
(9) 

(10) 



where in the first line, is the projection into fixed 

total number of particles, i.e., holons and Nj spinous; 
while is the projection into physical Hilbert space, 
i.e., removing all states not satisfying constraint Eq.(3). 
In the second line, P/v is the projection into fixed to- 
tal number of electrons, which has to be Nf, Po is the 
projection which removes all double occupancies. a{k) is 
defined as a{k) = —. 

Projected BSC wavefunction turned out to be a sur- 
prisingly good trial wavefunction for t-J model[17]. How- 
ever numerical studies[18, 19] showed that the quasi- 
particle weight is almost a constant along the direction 
from (tTjO) to (0, tt), i.e., it fails to reproduce the di- 
chotomy. The quasi-particle current of pBCSwf is also 
pretty smooth in the k space[18]. It is because pBCSwf 
is unable to capture the momentum dependence proper- 
ties that we need a new approach to underdoped high 
superconductors . 



B. How to capture fc-dependence features? 
— Spinon-dopon approach and projected 
spinon-dopon wavefunction 

Rebeiro and Wen [20] developed this new mean- field ap- 
proach trying to capture the spinon-holon recombination 
physics. In the following we briefly review their work. We 
know that at low temperature, spinon and holon recom- 
bine pretty strongly to give electron-like quasi-particle. 
So it is natural to introduce dopon operator - a bound 
state between a spinon and a holon - to describe low 
energy excitations. Note that a dopon has the same 
quantum number as an electron and describes a doped 
electron (or hole). But the Mott and spin- liquid physics 
at half filling should also be addressed. So one should 
also keep the spinon operator. As a result, two types 
of fermions are introduced here: spinon / and dopon d. 
Spinon carries spin 1/2 and no charge, and dopon carries 
spin 1/2 and charge 1. By introducing these two types 
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of fermions one enlarges the Hilbcrt space: now there arc 
16 states per site, among them only three are physical. 
The three physical states on site i can be represented in 
terms of spinon and dopon fermions as: 

IT) = IT/). U) = U/). m = ^\hU-lfU). 

(11) 

Here please notice that the constraints are two-fold: 
firstly there must be one / spinon per site, secondly d 
dopon has to form a local singlet with the spinon. 

One can do a self-consistent mean-field study. The 
mean-field Hamiltonian takes the form: 

Hmean = {-'^Xicosk^ + COS ky) - IJ.)fljka 

+ A(cosA;x - COS ky)fl^fi^^ + efe4a4c« 

+ PkfLdka + h.c.. (12) 

Here Hmean can be divided into three parts: spinon 
part, dopon part and spinon-dopon interaction. The 
spinon part describes the usual rf-wave paired ansatz: 

X = JifljkJ' ^ = J{fk]f-ki)- The dopon part is 
simply a free dopon band, with efc determined by high en- 
ergy ARPES measurement. Note that is not taken as 
tunable mean- field parameter. Finally the spinon-dopon 
interaction is described by a fc-dependent hybridization, 
roughly speaking /4 = ek{d\j^). One can see that d\f^ 
is a bosonic field carrying charge 1 and spin 0. Its non- 
zero average value corresponds to holon condensation in 
slave-boson approach, which leads to superconductivity. 
fi is the chemical potential required to time the doping. 

Along this line Rebeiro and Wen did a mean-field phase 
diagram, and successfully fit to ARPES data and tun- 
neling data [23, 24]. Here we try to emphasize that the 
main lesson we learned from this new mean-field ap- 
proach is that one can have a fc-dependent hybridization 
at mean-field level (in Eq.(12) this hybridization is con- 
trolled by f3k and energy spectrum of spinon band and 
dopon band.), which is roughly the coimterpart of holon 
condensation in slave-boson approach. This is why one 
can study detailed features like dichotomy in this new 
approach. 

Several open questions naturally arise in this new ap- 
proach. It seems there are two types of excitations, 
spinon and dopon, what do they look like? We also know 
that mean-field approach is not very reliable, so it would 
be nice to understand the physical trial wavefunction cor- 
responding to the new mean-field approach, from where 
we would know exactly what we are doing. In the follow- 
ing we try to answer these questions. 

Let us construct trial wavefunctions based on this 
spinon-dopon mean-field approach. One can simply take 
a mean-field ground state wavefunction, then do a pro- 
jection back into the physical Hilbert space, just like the 
way we did in the slave boson case: 

\^psd) = PsdPn\^sd .mean )■ (13) 



Here P/v is the projection into fixed number of spinon 
and dopon, which gives the correct doping; and Psd 
is the projection into to physical Hilbert space Eq.(ll). 

mean) IS the ground state wavefunction of some 
mean-field Hamiltonian in the form of Eq.(12). Sup- 
pose we know how to do this projection numerically, one 
can do a variational study of the these Projected Spinon- 
Dopon Wavefunctions (pSDwf), to see what is the lowest- 
energy ansatz. In general, however, the full projection is 
not doable, so we develop a simple numerical technique 
to do a local projection to have some rough idea about 
what kind of wavefunction is energetically favorable (See 
Appendix A). What we found is that the best trial wave- 
function for underdoped case has the following form: 

I^Isd) = PsDPN\^fh,mean) (14) 

= PsdPn exp b{k)jyiu^ |0) (15) 



DC 




where 

/L = yW^/L + /3;c4.- (17) 

Here / form a d-wave paired state and the superscript 
SC means this wavefunction is superconducting. h{k) 
and /3fe are some real functions and we assume /3fe = P-k 
to respect time reversal symmetry. For this particular 
ansatz, full projection is doable in low doping limit. In 
sciction IV wc develop the numerical method to do the 
full projection and we will see that this wavefunction is 
a even better trial wavefunction than pBCSwf. 

Note that the total number of / and d fermions is + 
Nh- Also Psd requires one /-fermion per site, so totally 
N /-fermions. Therefore we must have Nh rf-fermions, 
which gives the correct doping. 

C. How does pSDwf capture the fc-dependent 
features? — properties of wavefunction before 
projection \^sD,mean)'- at mean-field level. 

The form of |$P5^) looks very similar to pBCSwf, ba- 
sically we arc constructing a pairing wavefunction based 
on hybridized fermion fk ■ In the next section we will see 
that l^p^jj) and pBCSwf are indeed closely related. For 
the moment let us have a closer look at the wavefunction 
\^SD mean) bcforc projcction. The idea is that physical 
properties may not change drastically after the projec- 
tion. In this case the mean-field level understanding will 
give us insight of the wavefunction after the projection. 

First of all it is obvious that this wavefunction is su- 
perconducting. That is because the nonzero Pk signals 
the mixing between spinon and dopon (fl^dka) 0) and 
thus signals breaking of charge conservation. It is natural 
to believe the superconductivity survives after projection. 
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Let us introduce the other combination of / and d 
fermions; 

4a = -/3fe/L + ^1^4a> (18) 



and the quasi-particle operators: 



where 



Uk 



1-kl = Ukf-kl + Vkflp 

m 



1 



Vk = 



(19) 
(20) 

(21) 



are the coherent factors for a rf-wave paired state. We 
can show that 

= exp Ij2 mfl^fUi] |0), 



satisfies: 



and 



d'ka\^SD,mean) — 0) 



(22) 



lka\^SD.,mean)^0- (23) 

The mean-field Hamiltonian which can generate 
\^SD,mean) ^ ground State is simply: 

Hmean = fif^hka^-ka + ed(^)4«4a) - (24) 

k 

with e^{k),eg{k) > 0. Later we will see that there are 
physical reasons that ej(A;) > ef{k), meaning / band 

is lowest energy excitation, and d band is fully gapped, 
ej(fc) > for any k. 

We can express / and d fermions in terms of 7 and d 
fermions: 



fka = \J'^- 0k{^klka + Vk^aSl-u) " l^kdka (25) 

dka = PkiUklka + VkeaSliks) + y'' 1 " /3fe4a- (26) 

Based on Eq.(25) and (26), it is easy to obtain: 

(^SC Iff f |^5C \ 

\ SD.raean\J ka-> ka\ SD.nieanl /o2\~2 /orrN 

™ ^Tsc ^ = (1 - l^k)^k (27) 



\^ S D ,mean\\^ S D .mean/ 
^SD,mean\^ka"'ka\^SD,mean/ 

~1¥ 



sc w^sc \ 

SD,mean\\ SD^meanl 



PkVl 



(28) 



We know that the mean-field wavefunction should give 
one /-fermion and a; = ^ rf-fermion per site on average: 



k 



(29) 
(30) 



In the low doping limit a; ^ 0, it is clear from the above 
relations that 0^. oc x. 

Now let us understand how to calculate Z_ and Z+ on 
this mean-field wavefunction. and Z- are defined to 
be: 



Z-,k 



■ ■ 



\{N-l,k\ck\^^s 



(AT. 



l,k\N 

\{N + 1J,\( 



-GS\^GSI 



{N + l,k\N + l,k){^^s\ns) 



(31) 

(32) 



where |7V — 1, fc) (|Ar-|- 1, A:)) are the lowest-energy iV — 1 
{N + 1) electron states which have nonzero overlap with 

In our mean-field wavefunction, the lowerest energy ex- 
cited states are given by creating 7fe-quasi-particlc. Note 
that now d}^ is the hole creation operator, so at mean-field 
level the Zk for spinon-dopon wavefunction are: 



ySD _ 



«l7.T7^Tl*fg,n.ean>(*ig,n.eanl*|g.„»eaj' 

= Pl4 (33) 

+ <*fg,„.eanl7_,i7l.J*ig,^ean>(*fg,„.eanl*|g.„.ea„> 



(34) 



At this moment, let us compare spinon-dopon wave- 
function (SDwf) with BCS wavefunction (BCSwf), both 
before projection (In Section II D we will compare them 
after projection). 

In Section II A we view pBCSwf as the projected slave- 
boson mean-field state into physical Hilbert space. We 
may also view pBCSwf as projected BCSwf with all dou- 
ble occupancies removed: 

l^PBCs) = PdPnI^bcs) (35) 
= PdPn exp(^ a{k)ci^cl^^)\0) (36) 



Nj/2 



|0), (37) 



where 



|$BCs) =exp(^a(fc)4cL,^)|0) (38) 
fe 

ocJliuk + Vkcl^cl^^m. (39) 
fe 

Before the projection, the spectral weight of the electron 
operator Ck can be calculated easily: 



yBCS _ .,1 _ „ 
^-,fe -^k-^k 



Zl%^ = ul = X- 



nk, 



(40) 



where = {c\ci.)bcs- For a rf-wave BCSwf Eq.(6), we 

can plot the Z^ in Fig.l. In low doping limit, parameters 
are taken as /x = 0, x = 1, A = 0.55. Such choice of 
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FIG. 1: Plot of (left) and (right) within one quarter 
Brillouin Zone, and ky range from to tt. 

parameters leads to a pBCSwf with lowest average energy 
at half filling. 

The Z+^k and for the pBCSwf after the projec- 

tion were also calculated [18] . Roughly speaking what was 
found is that the Z^. profile after projection is similar to 
that before the projection. There is a quasi fermi sur- 
face, which is roughly along the diagonal direction, Z-^k 
is large inside the fermi surface and decreases very fast 
when you go outside fermi surface; while Z+^k is large 
outside fermi surface, and decreases fast when you go 
into fermi surface. But there is one big difference, which 
is a reduction factor. For the Z^^k, this reduction factor 
was found to be proportional to x. But for Z_ k, this 
reduction factor depends on k and is finite (around 0.2) 
for k = (0, 0) even at low-doping limit. From slave-boson 
approach Eq.(4), we already see that Z (x x a,t mean field 
level. Basically at half filling, Z = and we have a Mott 
insulator instead of a band insulator. 

Notice that along diagonal direction (tt, 0) — (0, tt), the 
Z^'~^^ is dispersionless: Z^'~^^ = 0.5, which does not have 
dichotomy feature; After projection, there is a factor x 
reduction, but Z^ is still almost a constant along the 
diagonal direction[18, 19]. 

To compare the calculated Z from the BSCwf and 
SDwf, we note that the projected wave functions, pBC- 
Swf and pSDwf, are closely related (see section II D). 
More precisely: 

\^¥sd) ^ lapses) if: 
Uk = Vk, Vk = Uk, and Pk = Po (constant) (41) 

It is easy to understand this identification at half-filling, 
since both wavcfunctions simply give the same spin-liquid 
(usually referred to as staggered fiux spin liquid in liter- 
atures), characterized by x ^-^d A. Now in the pSDwf 
f = f, with no mixing with d-fermion. It is simply a 
particle-hole transformed pBCSwf, by which Uk trans- 
formed into Vk and vice versa. The important message 
is that Uk,Vk and Uk,Vk characterize the spin dynamics, 
but Pk characterizes the charge dynamics. 

With this identification in mind, from Eq. (33,34) and 
(40) we immediately know that when Pk = Po these two 



wavcfunctions give the same mean-field Zk profile except 
that SDwf has an extra x factor, because Pf — Pq oc x 
in low doping limit. However, when Pk has a strong k- 
dependence, Zk from the two approaches can be very 
different. 

Let us think about whether or not these wavcfunctions 
can capture dichotomy in low doping limit. What did we 
learn from these mean-field result? We learned that it 
is impossible to capture dichotomy by BCSwf, because 
in order to capture the fc-dependence along (0, tt)-{tt, 0), 
one has to tune m^, Vk- Because d-wave Uk,Vk are con- 
stant along (0, 7r)-(7r, 0), one has to destroy the d-wave 
ansatz to have a fc-dependent Uk, Vk along (0, 7r)~(7r, 0). 
This leads to a higher J energy. On the other hand, it 
is possible to capture dichotomy by SDwf, because one 
can tune Pk to have a strong fc-dependence while keeping 
Mfe, Vk to be d-wave ansatz. This will not destroy the spin 
background. Based on our experience of projection, we 
expect that even after projection, the above statement is 
qualitatively true. 



D. Why mixing Pk has a strong fc-dependence? 
—Relation between pBCSwf and pSDwf 

In the last section we see that pSDwf can potentially 
capture the dichotomy through a k dependent Pk- Now 
the issue is, why does the Pk want to have a k dependence 
that can explain the dichotomy in Z^? Why does such 
fc-dependent Pk lead to a pSDwf which is energetically 
more favorable? To understand this, we need to know 
what a pSDwf looks like in real space. 

The discussion below for identifying the relation be- 
tween pSDwf and pBCSwf is rather long. The result, 
however, is simple. Let us present the result here first. 
We introduce Pk — Pk/ \/l — Pi- In low doping limit, 
Pk ~ Pk- For the simplest one-hole case, if Pk has the sim- 
plest modulation in fc-space Pk = /3o-|-2/3i(cos fca;-|-cos ky), 
then the pSDwf can be viewed as pBCSwf mixed with 
the wavefunction generated by the nearest neighbor hop- 
ping operators (see Eq.(56)). For more complicated Pk — 
Po + 2Pi (cos k^ + cos ky) + 4:P2 cos k^ cos ky + 2P3 (cos 2fca; -|- 
cos 2fcj,), the pSDwf can be viewed as pBCSwf mixed with 
the wavefunction generated by the nearest neighbor, next 
nearest neighbor and third nearest neighbor hopping op- 
erators (see Eq.(59)). Therefore to lower the hopping 
energy, finite /J^'s are naturally developed. This is why 
Pk with a proper fc dependence is more energetically fa- 
vorable. 

Before we look into pSDwf, let us review what a pBC- 
Swf looks like in real space. One can do a Fourier trans- 
formation: 
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l^PBCs) = Pd [^a{k)clcl^A \0)=Pd[ J2 " R-»T)4„,T4,i, J |0). 



(42) 



where a(r) = J2k ^i^) cos(fc-r). If we have a spin basis {Rj|,Rjj^}, where Rj| labels the positions of spin up electrons 
and Rji labels the positions of spin down electrons: 



({RiT,R,J|$PBCs) = 



a(R-ii~R-iT) a(Ri4 - R21 



a(R2x— Ri|) a(R2x — R2t) ••• a(R2i — R«/^) 

2 I 



a(R^|-RiT) o(R-£L/ , - ^2t) ••• a(R«v -RiVf ) 

2 4- 2-1- 2-1-2 1 



(43) 



We see that the overlap between a spin basis and pBCSwf is simply a single slater determinant of a two-particle 
wavefunction. This is why pBCSwf can be numerically simulated on a fairly large lattice. 
Now we go back to pSDwf. Up to a normalization constant, one can express pSDwf as: 



|0) 



(44) 



where (ik = Pk/ — Pi- Since Pk oc \/x, in the low doping limit, Pk = Pk- 
One can also do a Fourier transformation into the real space: 



where Ps's are the Fourier components of Pk- 

Pk=Po + P^e'^" + P-^e-'^' + pye'^y + p-ye'"''^ +■■■ . 

We should only consider the rotation invariant Pk, and let us only keep the first three Fourier components: 

Pk =Po + 2/3i (cos kx + cos ky) + AP2 cos kx cos ky 
+ 2p^{cos2kx + cos2ky). 



|0), 
(45) 

(46) 



(47) 



We claimed that if Pk = Po, then pSDwf is identical to pBCSwf if b{k) = Let us see how that is true. Without 
/3i,2,3, Eq.(45) is: 



\^'p%{Po)) =Psd\ E - ^Mk,^ + ^o4nT)(4u + 



|0). 



(48) 



What does a pBCSwf look like? If one does a particle- hole transformation c-^ hii , pBCSwf is: 



pd I E K^Ji-^^^)hk,Aui 



(49) 

|0), 
(50) 
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where Pd is the projection forbidding any empty site. 

If we consider a spin basis {Ri|, R-jx}, with the empty 
sites {Rfe^oji then after particle-hole transformation, we 
have single occupied sites {Ri|,Rjx}, and double occu- 
pied sites {Rfe_o}- So the position of spin up and down 



sites in the hole representation are {Ri|, Rj|};i, where 

{fi^^}h = {R^T}U{Rfc,n} and {fijijh = {Rjj,} U {Rfc,o}. 
The overlap of pBCSwf and the spin basis in hole repre- 
sentation is: 



({R,l,R^J|4>PBCs) = ({R,T,Rji}h|$PBcs) 

6(Ru-Rit) 6(Ru-R2t) ••• 6(Ru-R«i^ 



&(R2i - RiT 



6(R2x - R2t) 



5(R2i - R 



(51) 



b(R ]v+jv^ — Rit) &(R 



R2t) 



6(R 



RjV + Wj, ) 




FIG. 2: pSDwf with only Po- The site is empty, /-fermions 
are represented by green spin, and d-fermion is represented by 
red spin. Black valence bonds are bonds with of /-fermions, 
while red valence bond has a d-fermion. The two figures 
are two contributions of the overlap between pSDwf and a 
spin basis |l|3|2j^4x0emp), and they correspond to the same 
term in determinant Eq.(52). The two figures give rise to 
states: left: | Ti/io/T3/i2/TodJ,4/> = | Ti/T3/i2/i4/ (J,o/Tod)> 
and right: | Ti/iodT3/i2/To/i4/> = ~\ Ti/T3/i2/i4/ (To/iod)> 
The minus sign means that the two figures contribute addi- 
tively. 



The equation works this way because if one simply ex- 



pands the polynomial in Eq.(50), each sum will give you 
one term in the expansion of the determinant in Eq.(51), 
and Pauli statistics is accounted by the sign in determi- 
nant expansion. 

Now we can do the same analysis on pSDwf Eq.(48). 
First of all Pq is not relevant in the wavefunction, since 
we are projecting into a state with fixed number of d- 
fermion, which means that all f3o does is to give an over- 
all factor /3(^'' in front of the wavefunction. To have 
an overlap with spin basis {R^i, Rj|, Rfe^o}; we know 
that on empty site Rfe^o the expansion of polynomial 
Eq.(48) should give either | t/id) or | Ulf)- After 
projection each case would contribute to -^10) where 

\o) = ^{\hid) + \Uif)) 

One immediately sees that the expansion of polyno- 
mial Eq.(48) gives similar terms as the expansion of 
Eq.(50); actually corresponding to one term in determi- 
nant Eq.(51), we have 2^'' terms from Eq.(48), since we 
can either have | t/id) or | tdi/) for each empty site. 
The details are visualized in Fig. 2. Taking into account 
the factor-^ of projection, one has: 



({R.T,R,JI<i>^§D(/3o)) = (^2/30)^" 



6(Ru - Rit) 6(Ru - R2t) 
&(R2i - Rit) fe(R2i - R2t) 



6(R 



Rit) 6(R 



R2t) 



6(Rii — R jv+JVh - 

^ -.2 

^(R2i — R w + Wh - 



b(R w + w,, — Rw + jVfc ) 

2 ^ 2 



, (52) 



we found that \^pbcs) and l^p^oiM) are the same wavefunction. 
Now let us put in the simplest fc-dependence in (3k'. 



= /3o + 2;0i(cos kx + cos ky), 



(53) 



9 



we try to write \^p^jj{Po, in real space. After the Fourier transformation into real space: 



-psd\ E HR.u-RiT)(4,T+^o4.,T+^i E 4.,+„T)(4u+^o4,ii+A E <,+..i) 

(54) 



If we expand this polynomial Eq.(54), of course we will still have contribution from Pq terms which is nothing but 
the right hand side of Eq.(52). But apart from that, we also have contribution from Pi, which makes the problem 

more complicated. To start, let us consider the case of a single hole Nh = 1. To have an overlap with spin basis 
{Ri|, Rj|, R/c,o}) the d-fermion on empty site Rfe,o can also come from a bond connecting a spinful site and R/c,o + 
which is the /3i term effect. Let us consider the case S = y. We can also assume the spin state on site Rfc,o + S is spin 
down. Now it appears that we have two ways to construct the empty site on R^.q: | T/id) or I Tdi/)- We study the 
two cases separately. Firstly if the empty site is constructed by | t/id), shown in Fig. 3, careful observation tells us 
that the contribution to the overlap is exactly cancelled by fermion statistics. On the other hand, if the empty site is 
constructed by | t^i f)^ we have the case in Fig. 4. After careful observation, we know that this type of contribution is 

— ■^ times the overlap between |$po£)(/3o)) and the spin basis that differs from {Rjt, R,i , R^ o} by a hopping along 
y. Considering the fact that the shift can also be —y and ±x, one has: 



( {Ril , RjI , R/c,o } I *PSD (/^o , /3i )) 
=\/2^o({RiT.R'ji.R'M}l*Pscs) H 



V2 



({Ri1-,Rjj^,R/j,o}| 



E 

<5— ±x,±y 



=Rfc,o+5,aCRfc,o,al^PSCs), 



(55) 



where the minus sign in the second terms comes from Fermi statistics. 
Just by looking at Eq.(55), we arrive at the conclusion: 



\^'p%0oJl)) = I^PBCS) + ( ) Pd E 4+6,aCi,a\^PBCs). 



(56) 



Let us study Eq.(55). With out (3i, one has a single Slater determinant for the overlap with a spin basis; with /3i, 
we have 1 + Ushi /t = 5 Slater determinants, where Ushi ft is the total number of ways that one hole can hop. Later we 
will see that for nh holes, the number of Slater determinants for the overlap is (l + rishift)"'* , which means numerically 
one can only do few holes. 

The result (56) is obtained by studying one-hole case, and it is not hard to generalize the result for the multi-hole 
case. Basically, each hole may either not hop or hop once with a prefactor -^^j but not hop more than once. For 
example, for two-hole case, we have: 



I*psd(/3o,/3i)) 



-Pi J 

2/3o 



E 



V 



2po 



S"+i52 ,a2 ^',"2 '-'i+Si ,ai'^i,ai 



1$ 



/ 



PBCSl 



(57) 



where the constraint i ^ j + 5^ makes sure no hole can hop twice, and the coefficient ^ comes from double counting. 
We can also easily generalize it to the case with $2 and P^.... For two hole case, we have: 



*psd(/3o,/35)) 



,t 



^4 E 



-Ps, 



V 



61,62 



2/3o y V 2/3o 



\^PBCS)- 

(58) 
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In the end, the general formula for multi-hole pSDwf is: 



l^lfsoi^oJs)) = Pd exp„,„^=o,i I 1 + E -^4+5.0.0^,0. I Wpbcs), 



(59) 



where Uhop = 0, 1 ensures that no hole can hop twice. 

















% 



FIG. 3: pSDwf with only The site is empty. /- 

fermions are represented by green spin, and d-fermion is rep- 
resented by red spin. Black valence bonds are bonds with 
of /-fermions, while red valence bond has a d-fermion. No- 
tice that the position of d-fermion is shifted by — y by /3i 
term effect (red dotted line). The two figures show the two 
contributions of the overlap between pSDwf and spin basis 



|lT3T2i4|0e 



with spin of f-fermion on site 2 and spin 



of d-fermion on site are parallel. They give rise to states: 
left: I Ti/iodT3/i2/To/i4/) = I Ti/T3/i2/i4/ (iodto/)) and 
right: | Ti/i2/T3/iodTo/i4/> = -| Ti/T3/i2/i4/ (iodto/)) 
The minus sign means the two figures contribute subtrac- 
tively, i.e., they cancel exactly. 



III. HOW TO MEASURE THE MIXING /3fc? 
-PHYSICAL MEANING OF SPINON 
EXCITATION AND DOPON EXCITATION 

From Eq. (34,33), we know what at the mean- field level, 
f3k can be measured by spectral weight Z . After projec- 
tion, it is natural to expect that Z is also closely related 
to Pk- But how to calculate Z after projection? Basically 
through Eq.(32,31), we need a good trial ground state 
and excited state. The ground state would be nothing 
but pSDwf. What is a good excited state? To be spe- 
cific, let us study Z_, then the question is how to obtain 
|iV- 1)? 

In pBCSwf, the good excited state (referred as quasi- 
particle state) is found to be: 

1^^ - l>f/'''^ = PdCpPn\BCS), (60) 

where Pn project into fixed Nf number of electrons. The 
way we construct excited state here is simple: first find 
a excited state on the mean- field level, then do a projec- 
tion. In pSDwf, which includes pBCSwf as a limit, we 
should have a similar formula. But now we have two pos- 
sible ways to construct excitation states, since on mean- 
field level we have two types of fermions / and d, they 
correspond to two types of excitations. Now it is impor- 
tant to understand what each type of excitations looks 







'l \ A 













FIG. 4: pSDwf with only Pi. The site is empty, /-fermions 
are represented by green spin, and d-fermion is represented by 
red spin. Black valence bonds are bonds with of /-fermions, 
while red valence bond has a d-fermion. Note that the po- 
sition of d-fermion is shifted by —y by /3i term effect (red 
dotted line). This figure shows another contribution to the 
overlap between pSDwf and spin basis |l|3|2|4|0em,p), with 
spin of f-fermion on site 2 and spin of d-fermion on site are 
anti-parallel. It gives rise to state: | Ti/io/T3/i2/Totii4/) = 
~l Ti/T3/i2/i4/ (Todio/))- Note that there is a contribution 
for (lT3T0i4 l^'psoiPo)), with the same valence bond 
map, where |l|3|0x4x2emp) is the result state after a hopping 
along y acting on original state |l|3|2j4|0emp). That one 
would give a state | Ti/io/T3/i2/T2di4/) = | Ti/T3/io/i4/ 
(T2di2/))- The minus sign in the shifted d-fermion overlap 
comes from fermi statistics. 



like. It turns out that the f-type excitation corresponds 
to the quasi-particle excitation, and d-type excitation cor- 
responds to bare hole excitation. Thus the quasi-particle 
state I A'' — 1) for calculating Z_ is the /-type excitation. 

The main result for this section is Eq.(71,72) for /- 
excitation and Eq.(77) for d-excitation. One can see that 
the /-excitation of pSDwf is just the quasi-particle exci- 
tation in pBCSwf together with hopping terms acting on 
it. And c?-excitation is the bare hole on a pSDwf ground 
state. Let's see how those happen: 

For /-excitation, 



\N-l)f = PsDflpPN\^l%) 



(61) 



Here \N — 1) f has Nf — 1 number of electrons, because 
before projection there are totally N -\- Nh -\- 1 fermions, 
and we know projection enforces one /-fermion per site, 
so totally there are N/^ -f 1 d-fermions, i.e., holes. What 
does this wavefunction look like after projection? 
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To understand this we first try to understand the ex- 
citation of pBCSwf. What is the excited state in terms 
of spin basis? One way to see it is to identiiy: 

\N - l)f/^^ = Pnc,^ (J2 «(^)4TcLfe J |0) (62) 

^^^^-Pi(E«w4TcUJ io). 

(63) 



In this form the overlap with a spin basis {Ri|,Rj|} is 
easy to sec. Notice now the number of up spins is Nf/2 — 
1, and the number of down spins is Nf/2, so there is one 
more site in {Rjj^}. The only fashion to construct a spin 

basis is: let clp^ create an electron somewhere, then let 

a(fc)cj[.^clj,| create the valence bonds. After observation, 
the overlap is: 



({R,T,R,J|iV-l)f/^^ 



"(R-U^R-it) «(R-ii-R-2T) 
a(R2i - Rit) a(R24 - R2t) 



a(Ri|-Rjv^_^P e" 
a{R2i - Ri_i|p e-'P^-i 



(64) 



a(Ri^i -R-it) a(R«, 

2 2 



R2t) 



ipH. J 



But to compare with pSDwf formalism, we want to see the same result in a different way. Let us do a particle-hole 
transformation, just like what we did in Eq.(50). 



Nf/2 



(N+Nh)/2 



\0) cx Pnhl^^ (J2b{k)hlh^ 

V k 



|0). (65) 



The only way to construct a spin basis in hole representation {R-if, Rjj.}/i is to let h^_p^ construct a hole somewhere, 
and let 6(A;)/i^^/it^j^ construct the valence bonds to fill the lattice. The overlap in hole representation is: 



({R,T,R,J|Ar- l)f/^^ = ({R,T,R,J,|iV- l)f/^^ 



6(Ru - Rit) 
6(R2i - Rit) 



6(Ru -R2t) 
6(R2i - R2t) 



6(Ri1 -Riv 



6(R2j. — Rn+i 



6(R iv+jv^ ^ - Rit) b(R iv+w^, , - R2t) • • • 



g-jpRij, 
g-ipR.2i 



(66) 



Now let us go back to pSDwf. What is an /-type 
excitation? The only way to construct a spin basis, 
|{Ri|, Rjt}^), is to let /Ip construct an /-fermion some- 
where, then let b{k)(fl^ + Pkd\^){f^_ui + hd^^^,^) con- 
struct the valence bonds to fill the whole lattice. We first 
consider the case $k = $0- In this case, the constructing 
precess is in exactly the same fashion as in Eq.(66), ex- 
cept for one difference: there is a coefficient of V2(3q for 
each hole. That is because for each hole there are two 



contributions, one from | ]fid), the other from | tdi/); 
each with coefficient of -\/l/2/3o; unless that the hole and 
the spinon created by /l^ are on the same site, in which 
case we have only one contribution. If we ignore the last 
effect (since it is an infinitesimal change to the wave- 
function in the low doping limit, and it also comes as 
an artifact of our projective construction), we conclude 
that: 
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({RiT,RjJ|iV- l,^o)/ = ({Rn,RiiU|A^- 1,^)/ 
6(Rii - Rit) b{Rii - R2t) 

6(R2i - Rit) 6(R2i - R2t) 



6(R2i - Rw+iv^ ) 

2 I 



6(R iv+iV;. ^ - Rit) ?>(R w+iVfe ii^j^ ~ R2t) ••• b(R w+iVh I ^^j^ - R w+iVfe P e 

PBCS 
qp ) 



g-ipR2i 



"'P"-'^+/fr+i,i 



SO 



|7V-l,/3o)/cx |iV-l) 



PBCS 



(67) 



(68) 



(69) 



(70) 



The point is that /-type excitation describes the spin-charge separation picture of the excitation, because the hole 
and the unpaired spinon created by /l^ can be arbitrarily separated. And it turns out to be the low energy excitation 
of t-J model. 

What if the mixing (3k has momentum dependence? Similar to our study for the ground state wavefunction leading 
to Eq.(55,56), one can convince oneself that, in the one hole case 



PBCS 
qp 



-P2 

Wo 



D 



i,5=±(x+y),±(x—y) 

And for multi-hole case, similar to Eq.(59) 



2/3o 



PBCS 
qp 



(71) 



i,5=±2x,±2y 



PBCS 
qp 



(72) 



For d-excitation, story is different. It turns out d-excitation corresponds to bare hole excitation. What is a bare 
hole excitation |A'' — l)bh! For a pBCSwf, 



\N -l)^^'''' = Cj,^\^PBCs), 



(73) 



in terms of spin basis, it is easy to show that: 

{{Iii^,-Rjl,-Rko}\N -1)^^^^ = {{Ri^,-Rj^,-Rko}\cp^\^PBCs) oc ^e-'f^'=«({{RiT,Rfeo},R, J|$PBcs) 



E- 



a(Ri|-Ri|) a(Ri|-R2t) ••• a(Ri| - Rnj. ) a(Rij, - Rfco) 

2 -'- 

a(R2i-RiT) a(R2i-R2T) ••• o(R-2i - R«/ J a(R2i - Rfeo) 



a(Rjv^ — Rif) a(Rjv^ — R2-f) ••• a(Rjv^ — Rjv^_ ) a(Rjv^ — R^^o) 
What is a d-type excitation in terms of spin basis? We first consider the case where jSk = Po^ 



(74) 



One can see that the only way to construct a spin basis, |{Ri|, Rj|}/(), is to let &(fc)(/^| ~^ (^kdi-^){f-ki + PkdLkO 
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construct the valence bonds to fill the whole lattice, then 
find a site occupied by one /|-fermion only, and let dLp^ 
construct a hole there. By observation, we conclude, 

|7V-l,/3o)d(x|iV-l)f,^^^. (76) 



If Pk has /c-dependence, one can also convince oneself 
that 

\N-lJo,Ps)d (X Cp\<^lfsD{Po,m = \N-1)^^^^, 

(77) 

i.e., d-type excitation corresponds to the bare hole in 
pSDwf. 

To summarize, we have the following identification: 
/-type excitation corresponds to the low energy quasi- 
particle excitation, i.e., a state constructed by putting 
Cp operator inside the projection; d-type excitation cor- 
responds to the bare hole excitation, i.e., a state con- 
structed by putting Cp operator outside the projection. 



IV. NUMERICAL METHODS AND RESULTS 

We use Variational Monte Carlo (VMC) method to 
calculate the ground state energy (of 2 holes), the excited 
state energy (of 1 hole) of pSDwf and pBCSwf and the 
spectral weight Z-^k. 



Our pBCSwf calculation is mostly traditional. Never- 
theless the previous calculation of Z_ [19] is indirect and 
having uncontrolled error bars inside the fermi surface. 
We developed a straightforward technique to calculate 
Let us recall the definition of Z- Eq.(31). For pBC- 
Swf, if we relabel |iV - l)^,^^^ as \bh) and \N - l)^^^^ 
as \qp) to save notation: 



\{qp\hh)\-' 



{qp\qp){'^PBcs\^PBcs) 

\{qp\bh)\'' {bh\bh) 
{qp\qp){bh\bh) {^pbcs\^pbcs) 

\{qp\bh)\^ 



{qp\qp){bh\bh) 



rik, 



(78) 



where is the occupation number of particles at mo- 
mentum k. Uk can be calculated by VMC approach 
pretty straightforwardly [25]. In particular, one can easily 
show that at low doping limit, which is the case consid- 
ered in this paper, = ^ independent of k exactly. The 
only thing one needs to worry about is the overlap pref- 
actor between \qp) and \bh). Instead of calculating the 
factor itself, one can split the calculation into two. If we 
denote a spin basis as |s). 



{qp\bh) 



= E 

S 

{bh\bh} = ^ 



{qp\qp) 

{bh\qp) 



{qp\s){s\bh) 

{qp\qp) 

{bh\s){s\qp) 
{bh\bh) 



= E 

s 

= E 



{s\bh)\{qp\s)\'' 
{s\qp) {qp\qp) 

{s\qp) \{hh\s)\' 
{s\bh) {bh\bh) 



E 

s 

E 



{s\bh) 
{s\qp) 

{s\qp) 
{s\bh) 



Pbh{s). 



(79) 
(80) 



Since both {s\qp) and {s\bh) are Slater determinant or 
sum of Slater determinants (see Eq.(64) and Eq.(74)), 
the above two quantities can be calculated by Metropolis 
program in a straightforward fashion. Then the product 
of the two gives the Z-^k- This algorithm works for finite 
doping case, too. 

For pSDwf, because we include /c-dependent mixing, 
in each step of Metropolis random walk, we need to keep 
track of all the (1 -|- nsUft)"'^ matrices, which limit the 
calculation for few holes. 



A. Ground state at half filling and 2 holes 

The calculation is done for t-t'-t"-J model on 10 by 10 
lattice, where t = l,t' = -0.3, t" = 0.2 and J = 0.3. We 

choose periodic boundary condition in x-direction, and 
anti-periodic boundary condition in y-direction. 



For variational parameters, we choose the lowest- 
energy ansatz in Eq.(7)[26] with parameters 

X = l, A = 0.55, M = 0. (81) 

The energy for half-filling ground state is listed in Table 
I. 

For two holes, we compare the energy of ground states 
of pBCSwf and pSDwf. For pSDwf, to lower the t hop- 
ping energy, since t <0,hy Eq.(59), the sign of /3i should 
be negative. Similarly since t' > 0, t" < 0, the parame- 
ters lowering t' and t" hopping energy have signs /32 > 
and /33 < 0. We did a variational search for the optimal 
values of The results are listed in Table II, where we 
also compare it with pBCSwf with longer range hoppings 
(see Section IV E). 

We find that the energy of the best pSDwf is lower 
than the energy of the best pBCSwf. We note that the 
pSDwf and pBCSwf are identical at half filling. So the 
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energy per bond 


Si ■ Si+1 


-0.1710±0.0001 


-0.3200±0.0001 



TABLE I: half-filling ground state on 10 by 10 lattice 



wavefunction 


A 
X 


xl 

X 


x" 

X 


/3l 
So 


/32 

00 


la 

00 


total energy 

per bond 


{Si ■ Si+i) 
per bond 


Ti 


T2 


n 


pBCSwf 


0.55 

















-0.1872±0.0001 


-0.2977±0.0002 


2.64±0.01 


0.52±0.01 


0.48±0.01 


pBCSwf(optimal) 


0.55 


-0.4 


0.0 











-0.1890±0.0001 


-0.2947±0.0002 


2.66±0.01 


0.06±0.01 


1.07±0.01 


pBCSwf 


0.55 


-0.5 


0.1 











-0.1885±0.0001 


-0.2872±0.0002 


2.66±0.01 


-0.23±0.01 


1.52±0.01 


pSDwf(optimal) 


0.55 








-0.3 


0.3 


-0.1 


-0.1918±0.0001 


-0.2943±0.0002 


2.86±0.01 


-0.46±0.01 


0.77±0.01 



TABLE 11: Two holes on 10 by 10 lattice, t = 1, t' ^ -0.3, t" = 0.2 and J = 0.3. Ti, T2 and T3 stand for nearest neighbor 
hopping per hole ■^J2i,s=±x,±y{4ci+s), next nearest neighbor hopping per hole ■w;;;J2i,s=±(S:+y),±(S:-y){cUi+s) and third 
nearest neighbor hopping per hole ^=±24 ±2y('^i'^i+s) respectively. We compare pBCSwf of d-wave ansatz, pBCSwf with 

longer range hoppings x' and x"j and pSDwf with non-local mixings. The best trial pSDwf has an energy 1.5% below that of 
the best trial pBCSwf with longer range hoppings. Comparing the first and the last line which have the same spin correlation, 
we find that the energy of a hole in pSDwf is 0.46t lower than that of a hole in pBCSwf. Note that pBCSwf with longer range 
hoppings destroys the d-wave spin background. As a result, the attempt to lower the hopping energy by tuning x' and x" is 
not effective since it would increase the spin energy a lot. 



energy difference between the two states is purely a dop- 
ing effect. Comparing the first and the last hnc in table 
II which have the same spin correlation, we see that the 
total energies of the two states differ by 0.0046 x 200 since 
the the 10 by 10 lattice has 200 links. This energy differ- 
ence is due to the presence of two holes. So the energy 
of a hole in pSDwf is 0A6t lower than that of a hole in 
pBCSwf. This energy difference is big, indicating that 
the charge-spin correlation is much better described by 
pSDwf than pBCSwf. 



B. Hole doped case, quasi-particle excitations and 

Z-. 

In this section we study the excitations of t-i!-t"-J 

model, which is one hole on 10x10 lattice. We also 
compare pSDwf with pBCSwf. We know from Eq.(72) 
that the pSDwf /-excitation state goes back to pBCSwf 
quasi-particle excitation state when all non-local mixings 
(is = 0. Also from Eq.(72), one can see that to lower the 
t, fl and t" hopping energy, we should also have /3i < 0, 
/J2 > and /J3 < 0. Actually in the low doping limit, one 
should expect the non-local mixing (5s for quasi-particle 
excited states (/-excitation) to be same as the ground 
state. Here we adopt the values of /J^ from our study of 
2-hole system ground state. 

Our VMC calculation shows that the pSDwf or pBC- 
Swf has finite Z- deep inside the fermi surface even in the 
low doping limit x — > 0. This is physically wrong because 
deep inside fermi surface there is no well-defined quasi- 
particle, and the idea of calculating Z_ by a single parti- 
cle excited state is also incorrect. Nevertheless, because 
the low energy excitation is more and more quasi-particle 
like as one approaches the fermi surface, we expect that 




FIG. 5: Shapes of functions cos k^c cos ky (up) and cos 2fca; + 
cos2/Cj, (down) along diagonal direction from (7r,0) to (0,7r). 



the Z_ calculation remains valid close to fermi surface, 
roughly speaking, along the diagonal direction from (tt, 0) 
to (0,7r). 

From Eq.(33), we know that at the mean-field level, 
the modulation of Z_ is controlled by 0\. It is im- 
portant to study the shapes of (5\ for various cases. 
In Fig. 5 we plot the shapes of functions cos fc^ cos fcy 
and cos 2fc;j; -|- COS 2fcy along the diagonal direction. If 
/3fe — /3o + 2/3i(cosA:j; + cosA:j,), remains constant along 
the diagonal direction. If (3). = Pq + 4/^2 cos cos ky, for 
small /32 > 0, is reduced at the anti-nodal point. If 
/3fe = /3o + 2/33 (cos 2A;j; + cos2ky), for small /S^ < 0, (Sf is 
enhanced at the nodal point and suppressed at the anti- 



15 



(0.7t) 



r 



(jt.ji) 



(0.71) 



0.116 0.090 0.035 0.006 0.000 0.000 0.027 0.028 0.020 0.005 0.000 0.000 

0.138 0.117 0.042 0.001 0.001 0.002 0.060 0.059 0.030 0.002 0.001 0.002 
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FIG. 6: For one hole on 10 by 10 lattice, we plot Z-^^ of 
pBCSwf (left, X = 1, A = 0.55) and pSDwf (right, x = 1, 

A = 0.55, = -0.3, |i = 0.3, = -0.1). pBCSwf has 

PO PO PO 

almost constant Z- along diagonal direction from (tt, 0) to 
(0, 7r); while pSDwf has Z- suppressed at anti-nodal point. 



nodal point. Let us remember this trend: positive (32 
and negative drive the modulation of Z_ in the way 
consistent with dichotomy for hole doped samples. 

For small values of f3k we know that f3k f3k- Eq.(33) 
suggests that along the diagonal direction 

Z- k oc {Po + 4/32 cos kx cos ky + 2/?3(cos 2kx + cos 2ky))^. 

(82) 



into hole-doped t-J model with parameters t, —t', —t", J, 
together with a (tt, tt) shift in momentum space. 

The approach outlined in Eq.(79) and Eq.(80) still ap- 
plies here. But because of the particle-hole transforma- 
tion, we are calculating Z+ of the original electron-doped 
system. Because t' > and t" < 0, to favor longer range 
hoppings, we must have $2 < and f3^ > 0, which dif- 
fer from hole-doped case by a sign flip. As a result, the 
Z-i- now will be suppressed at nodal point, but enhanced 
at the anti-nodal point. This is exactly what people ob- 
served in exact diagonalization[8]. 

We did a variational search for the optimal variational 
parameters for pBCSwf with longer range hoppings x' 
and x") and pSDwf with non-local mixings. In Table 
III wc compare the energy of pBCSwf and pSDwf with 
2 electron doped on 10 by 10 lattice. In Fig. 8 we plot 
the Z+ map of pSDwf, one can see pSDwf has spectral 
weight of anti-dichotomy shape. 

In Fig. 9 we compare the energy dispersion of one- 
electron quasi-particle excitations of pBCSwf and pSDwf. 
The energy of a doped electron in pSDwf is 0.25t lower 
than that of an electron in pBCSwf. 



D. A prediction 



But as a mean- field result, one should expect that the 
above equation is only valid qualitatively. In fact to 
crudely fit the relation of the modulation of Z- and (}, 
we found it is better to have some order of unity extra 
factor in front of (3 terms, and (3i also contributes to the 
modulation of Z- as a uniform shift. 



(X 00 + ^ + P2 cos kx cos ky + ^ (cos 2k., 



For t-J model without f/ and t" , there is no reason 
to develop a finite value of (32 and (33 since there is no 
longer range hoppings. As a result, one expects that Z_ 
remains almost constant along the diagonal direction. 

For J-t'-t"- J model with t>0,t' <0, t" > 0, we know 
that (32 > and (3^ < have to be developed to favor 
longer range hoppings. So one expect Z^ should develop 
dichotomy shape along the diagonal direction. 

In Fig. 6 we compare the Z_ ^ of pBCSwf and pSDwf. 
One can sec that pSDwf shows strong dichotomy. 

In Fig. 7 we compare the energy dispersion of one- hole 
quasi-particle excitations of pBCSwf and pSDwf. The 
energy of a doped hole in pSDwf is 0.38t lower than that 
of a hole in pBCSwf. 



C. Electron doped case 

In electron-doped case, one can do a particle-hole 
transformation, then multiply a (—1) for the odd lattice 

electron operators. By doing so, the original electron- 
doped t-J model with parameters t, t', t" , J transformed 



In hole-doped and electron-doped case, t > Q and t' 
and tf' have opposite signs, and as a result Z_ develops 
strong k dependence along diagonal direction. What if t' 
and t" have the same sign? If both t' > Q and t" > 0, 
one expects that (32 < and /^a < to favor longer 
range hoppings. But they drive the modulation of Z- in 
opposite ways. As a result, one expects that for certain 
la; -|- cos 2fcy))^. ratio of values oi t' > and t" > of order 1, their 
/o^\ effects cancel and Z- remains constant along the diagonal 
direction, but with an enhanced value of Z- than the case 
of pure t-J model. Similarly for certain ratio of values of 
t' < and t" < of order 1, Z- remains constant along 
the diagonal direction, but with a suppressed value of Z_ 
than the case of pure t-J model. These predictions can 
be checked by exact diagonalization. 



E. pBCSwf with longer range hoppings 

One can view pSDwf as an improved pBCSwf. We 
choose the d-wave pairing wavefunction b{k) with only 
nearest hopping x and pairing A parameters. Then (32 
and (3^ encode some second-neighbor and third-neighbor 
correlations. The price to pay is to include more than 
one Slater determinants in spin basis. One may natu- 
rally ask, suppose we insist working on pBCSwf, if one 
puts in longer range hopping parameters like x' and x" 
in the pairing wavefunction b(k), one also encodes some 
second-neighbor and third-neighbor correlations, which 
may lower the second-neighbor and third-neighbor hop- 
ping energies. But in this way one can still work with 
a single Slater determinant. If our pSDwf with no-local 




(0.0) 



(it.O) 



(7t,JC) 



(0,0) 



-36.5' 

(n.O) 



(II/2.7I/2) 



(0,1) 



FIG. 7: Quasi-particle spectrum for one hole on 10 by 10 lattice, t — 1, t' = —0.3, t" = 0.2 and J — 0.3. The black square 
shows the spectrum of d-wave pBCSwf with x ~ ^ S'Hd A = 0.55, and the red diamond shows the spectrum of pSDwf with 
X = 1, A = 0.55, = -0.3, 1^ = 0.3 and |^ = -0.1. One can see the first hole doped to {n/2,n/2). The energy of a doped 
hole in pSDwf is 0.38t lower than that of a hole in pBCSwf. 



wavefunction 


A 

X 


x' 

X 


x" 

X 


/3l 
/3o 


/32 
/3o 


Po 


total energy 
per bond 


{Si ■ Si+i) 
per bond 


Ti 


T2 




pBCSwf 


0.55 

















-0.1884±0.0001 


-0.2977±0.0002 


2.64±0.01 


0.52±0.01 


0.48±0.01 


pBCSwf (optimal) 


0.55 


0.2 


0.0 











-0.1888±0.0001 


-0.2964±0.0002 


2.61±0.01 


0.70±0.02 


0.20±0.02 


pSDwf (optimal) 


0.55 








-0.5 


-0.3 


0.3 


-0.1910±0.0001 


-0.2971±0.0002 


2.57±0.01 


0.86±0.02 


-0.72±0.02 



TABLE III: Two electrons on 10 by 10 lattice, t = 1, t' — —0.3, t" = 0.2 and J = 0.3, and we mapped it into a hole-doped 
model with t = 1, t' = 0.3, t" = —0.2 and J = 0.3 with a (7r,7r) shift in momentum space. Ti, T2 and Ta stand for nearest 
neighbor hopping per hole Ei,fc±4,±a(clci+i), next nearest neighbor hopping per hole Ei,i=±(4+ji),±(4-y) ) and 
third nearest neighbor hopping per hole Ei i=±24 ±2y{^\''i+^) respectively. We compare pBCSwf of d-wave ansatz, pBCSwf 
with longer range hoppings x' and x", and pSDwf with non-local mixings. The best trial pSDwf has energy lowered by 1.2% 
from the best trial pBCSwf with longer range hoppings. And comparing the first line and the last line which have the same 
spin correlations, we find that the energy of a doped electron in pSDwf is Q.2Qt lower than that of an electron in pBCSwf. 
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FIG. 8: For one electron on 10 by 10 lattice, we plot Z+^k 
of pBCSwf(left, X = 1, A = 0.55) and pSDwf(right, x = 1, 
A = 0.55, 1^ = -0.5, 1^ = -0.3, 1^ = 0.3). By particle-hole 
symmetry, the Z+ of one electron pBCSwf is identical to the 
Z- of one hole pBCSwf together with a (vr, tt) momentum 
shift. pBCSwf has almost constant Z+ along the direction 
from (tt, 0) to (0, tt); while pSDwf has Z+ suppressed at nodal 
point and enhanced at anti- nodal point. 



mixing is physically similar to pBCSwf with longer range 
hoppings, why should one bother to work with many 



Slater determinants? 

We want to emphasize that our pSDwf is physically dif- 
ferent from pBCSwf even after we include longer range 
hoppings x' and x" ■ We note that, in the infinite-lattice 
limit with a few holes, the pBCSwf cannot have longer 
range hoppings {i.e. x' — x" — 0). Otherwise we are 
considering some other spin wavefunction instead of d- 
wave wavefunction, which will increase the spin energy 
by a finite amount per site. Therefore x' and x" have 
to vanish in low doping limit. In contrast, for our pS- 
Dwf, the the spin energy is not affected by finite Pi in 
the zero doping limit. Thus in the low doping limit, the 
spin energy is perturbed only slightly by a finite (32 and 
(3^. On the other hand a finite P2 and (3^ make the hop- 
ping energy much larger than that of pBCSwf. So in the 
infinite-lattice limit with a few holes, Pi will be finite and 
the energy of one hole will be lowered by a finite amount 
by turning on a finite j3i. 

Physically this means that (32 and (3^ characterize the 
charge correlations, while x' and x" characterize the spin 
correlations. The above claim is supported by 2-hole sys- 
tem on larger lattice, i.e., by lower the doping. In Table 
IV we list the energies of pBCSwf with longer range hop- 
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FIG. 9: Quasi-particle spectrum for one electron on 10 by 10 lattice, t = 1, t' = —0.3, t" = 0.2 and J = 0.3 (one can map 
it into a hole-doped model with t — 1, t' = 0.3, t" = —0.2 and J — 0.3 with a (7r,7r) shift in momentum space.). The black 
square shows the spectrum of d-wave pBCSwf with x = 1 s^nd A = 0.55, and the red diamond shows the spectrum of pSDwf 



with X = 1, A = 0.55, ^ = 



doped electron in pSDwf is 0.25t lower than that of an electron in pBCSwf. 



-0.3 and ^ = 0.3. One can see the first electron doped to (7r,0). The energy of a 
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FIG. 10: For one hole on 10 by 10 lattice, we plot of 

pBCSwf with longer range hopping x' = —0.4 (x = 1, A = 
0.55). 



ping and pSDwf on 14 by 14 lattice. Comparing with 
Table II one can see the spin energy for pSDwf is lowered 
further than that for pBCSwf. 

Another way to see that these two wavefunctions are 
different is by calculating Z_. Numerical results show 
that pSDwf has dichotomy whereas pBCSwf does not. 
Actually on the mean-field level, a negative x' and/or 
positive x" even make the Z_ larger on the anti-nodal 
point than on the nodal point. After projection, we ob- 
serve that Z_ still remains almost constant along the di- 
agonal direction for pBCSwf with longer range hoppings. 
In Fig. 10 we plot the Z_ map of pBCSwf with longer 
range hopping x' = —0.4. 



V. CONCLUSION 

In this paper we studied a new type of variational 
wavefunction, pSDwf. It can be viewed as an improved 
pBCSwf, and the improvement is that pSDwf correctly 
characterizes the charge dynamics and the correlation be- 



tween the doped holes/electrons and the nearby spins. 
This physics was missed by the previous pBCSwf. As 
a result, pSDwf correctly reproduces the dichotomy of 
hole-doped and electron-doped Mott insulator. 

In pSDwf, we introduced two types of fermions, spinon 
/ and dopon d. Spinous / carry spin but no charge. They 
form a c?-wave paired state that describes the spin liquid 
background. Dopons d carry both spin and charge and 
correspond to a bare doped hole. The mixing between 
spinous and dopons described by /3o, 1^2 and /Sa leads 
to a d-wave superconducting state. The charge dynam- 
ics (such as electron spectral function) is determined by 
those mixings. f3o is the on-site mixing (or local mixing), 
and Pi, (32 and /?3 are non-local mixings corresponding 
to mixing with first, second and third neighbors respec- 
tively. If pSDwf has only local mixing, it is identical to 
pBCSwf. With non-local mixings, pSDwf corresponds to 
pBCSwf with hopping terms acting on it. Therefore the 
wavefunction develops finite non-local mixings to lower 
the hopping energies. In particular, for the hole-doped 
case, to lower t' and t" energies, the mixing is described 
by /32 > and Ps < 0. 

The pSDwf can also be obtained by projecting the 
spinon-dopon mean-field wavefunction into the physical 
subspace. Therefore, one expects that some properties 
of pSDwf can be understood from the mean-field theory. 
In the mean-field theory, it is clear that the modulation 
of in k space is controlled by the non-local mix- 

ings. Our numerical calculation of Z^^k shows that the 
above mean-field result is valid even for the projected 
wave function. We find that /?2 > and /Sg < give ex- 
actly the dichotomy of Z__fc observed in the hole doped 
samples. Because /?2 > and Ps < are driven by t' 
and t" , the dichotomy is also driven by t' and t" . Thus 
to lower the hopping energy, the spectral weight is sup- 
pressed in some region in fc-space. This result confiicts a 
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03 
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total energy 
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per bond 
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pBCSwf 


0.55 

















-0.1793±0.0001 


-0.3075±0.0002 


2.65±0.02 


0.43±0.02 


0.65±0.02 


pBCSwf 


0.55 


-0.4 














-0.1796±0.0001 


-0.3043±0.0002 


2.66±0.02 


0.02±0.02 


1.21±0.02 


pSDwf 


0.55 








-0.3 


0.3 


-0.1 


-0.1815±0.0001 


-0.3058±0.0002 


2.86±0.02 


-0.49±0.02 


0.87±0.02 



TABLE IV: Two holes on 14 by 14 lattice, t = 1, t' = -0.3, t" = 0.2 and J = 0.3. Ti, T2 and Ta stand for nearest neighbor 
hopping, next nearest neighbor hopping and third nearest neighbor hopping respectively. Although the spin energy of pBCSwf 
with finite longer range hoppings x' = —0.4 is slightly lower than that of pSDwf on 10 by 10 lattice, it is much higher on 14 by 
14 lattice. 



naive guess: to lower the hopping energy, the excitation 
should be more quasi-particle like. We also predict that 
the dichotomy will go away if t' and t" have the same 
sign and similar magnitude. In summary, we found a 
mean-field theory and the associated trial wavefunction 
capturing the dichotomy physics. 

Traditionally, in projected wavefunction variational 
approach, for example pBCSwf, people use wavefunc- 
tions which in real space correspond to a single Slater 
determinant. The reason to do so is simply to make the 
computation easier. Our study shows what kinds of im- 
portant physics that may be missed by doing so. In real 
space, the pSDwf is sum of (1 -|- rishift)"'''"''' number of 
Slater determinants, because each hole can either do not 
hop, or hop into one of Ushiji sites. So our calculation is 
limited to few-hole cases. However, the idea of introduc- 
ing many Slater determinant is quite general. For exam- 
ple, one can study another improved pBCSwf, which al- 
lows each hole to hop once but forbids two holes hopping 
together, therefore the number of Slater determinant is 
(1 -|- nshift)nhoie and many-hole cases are computation- 
ally achievable. This new improved pBCSwf is the first 
order approximation of pSDwf and remains to be stud- 
ied. For a long time there is a puzzle that doped Mott- 
insulator {ie the spin disordered metallic state) seems to 
be energetically favorable only at high doping x > 0.3. 
For X < 0.3 the doped spin density wave state have a 
lower energy. Our pSDwf may push this limit down to 
low doping which agrees with experiments better. This 
is because that including many Slater determinants can 
lower the energy per hole by a significant amount (about 
OAt). 

As we have stressed, pSDwf provides a better descrip- 
tion of spin-charge correlation, or more precisely, the spin 
configuration near a doped hole. This allows us to re- 
produce the dichotomy in quasiparticle spectral weights 
observed in experiments. The next question is whether 
the better understanding of the spin-charge correlation 
can lead to new experimental predictions. In the follow- 
ing, we will describe one such prediction in quasi-particle 
current distribution. 

We know that a finite supercurrent J, shifts the su- 
perconducting quasiparticle dispersion Efs- To the linear 
order in Jg, we have 

Ek{A)=Ek{0)+c-^jk-A, 



where c is the speed of light and we have introduced the 
vector potential A to represent the supercurrent: Jg = 
^^^^ A. jk is a very important function that characterizes 
how excited quasiparticles affect superfluid density ps- 
We call jk quasiparticle current. According to the BCS 
theory 

Jk = e-g^ = eWnormal, (84) 

where is the normal state dispersion which is roughly 
given by — —2t[cos{kx) + cos{ky)]. 

The previous study[18] of quasi-particle current for 
pBCSwf shows that the quasi-particle current is roughly 
given by the BCS result (84) scaled down by a factor a. 
Such a quasi-particle current has a smooth distribution 
in fc-space. Here we would like to stress that since the 
charge dynamics is not capture well by the pBCSwf, the 
above result from pBCSwf may not be reliable. We ex- 
pect that the quasi-particle current of pSDwf should has 
a strong fc-dependence, ie a large quasi-particle current 
near the nodal point where Zk is large and small quasi- 
particle current near the anti-nodal point where is 
small. Such a quasi-particle current distribution may ex- 
plain the temperature dependence of superfluid density 
[21]. 

Indeed, the mean-field spinon-dopon approach does 
give rise to a very different quasi-particle current dis- 
tribution which roughly follows Zk- For more detailed 
study in this direction and possible experimental tests, 
seeRef. [27]. 

We would like to thank Cody Nave, Sung-Sik Lee, P. 
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APPENDIX A: A SIMPLE ALGORITHM TO DO 
LOCAL PROJECTION 

Suppose the wavefunction before projection is the 
ground state of some fermionic quadratic Hamiltonian. 
One can always diagonalize the Hamiltonian so that all 
two-point correlation functions of fermion operators can 
be calculated exactly. For our SDwf, that means quanti- 
ties like (fifj), {d\fj),{d\dj)... can be calculated. 

Projection is supposed to remove the unphysical states. 
For a site i, the following operator removes the unphysical 
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states. 



Pi =nf,i{nf,i - 2){^n\i - ^nd,i + 1) 



(Al) 



It obviously ensures that njj = 1, Ud.i ^ 2 and / and d 
fermions form local singlet. To calculate energy, we do 
local projection on the relevant sites. For example, to 
calculate the J term energy, one actually calculates 



{Si ■ Sj)f 



{PiPj{Sf,i + Sd,i) ■ {Sf,j + Sd,i)PiP,) 



{Pipj) 



(A2) 



The denominator accounts for the wavefunction nor- 
malization due to projection. One can write operators 
PiSi ■ SjPj and Pi Si ■ SjPj in terms of fermion operators. 



By Wick's theorem, the expectation values of these op- 
erators reduce to a sum of products of fermion two-point 
correlation functions, which are known. Similarly for t 
term energy, one calculates for example. 



{'^i-\'^j'i)prj — 



{PiPjfl^h]hifj^PiPj) 
{PiPj) 



(A3) 



where h] = -^{fj^d]^ — fi^dl^) is the operator that cre- 
ates a hole at site i. 

One may ask whether we can do local projections on 
more and more sites, then the result will be closer and 
closer to the one of full projection. Unfortunately this 
cannot be done, because the number of terms in the sum- 
mation when we expand PiiPi2 ■ ■ ■ Pin increases exponen- 
tially fast as we increase n. Therefore we are limited to 
few sites. The above method can only be viewed as some 
renormalized mean-field approach. 
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